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III. REPRESENTATIONS OF PHOTON STATES 



1. Fock or "Number" States: 11 



As we have seen, the Fock or number states 



M) aII k.) 



[in-1] 



are complete set eigenstates of an important group of commuting observables — viz- 

Reprise of Characteristics and Properties of Fock States: 

a. The expectation value of the number operator and the fractional 
uncertainty associated with a single Fock state: 

(n\9t \n)= n [ III-2a ] 

An = [ "uncertainty" ] = ^{(n\^C 2 \n)-(n\9£ I")'} = [ m " 2b 1 

b. Expectation value of the fields associated with a single mode: 

For one mode Equations [ II-24a ] and [ II-24b ] reduce to 

E(? , t) = i e £ ja exp[/ k ■ r- /co f] - a ^ exp [-/ k f- i(0 t 



[ III-3a ] 



H(r, t) = i £ [k x e ja (t) exp [«' k • r- i to t - a \t) exp [-i* k • r- i to t j [ III-3b ] 



11 In what follows, for simplicity we drop the k subscripts on the operators and state vectors with the obvious 
meaning that |{n - }^ => | n) , a - => a , etc... 
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where £ = 




[n\E\n}= 
(w|H|w)= 



[ III-4a ] 



AE= l{(n\E-E\n)-(n\E\n) 2 } = ^(n + ±f = Jl t(n + ±f 





c. Phase of field associated with single mode: 

To obtain something analogous to the classical theory we would like to separate the 
creation and destruction operators (and, thus, the electric and magnetic field 
operators) into a product of amplitude and phase operators. Following Susskind 
and Glogower, 12 we define a phase operator, 4> such that 



Defined in this way, the basic properties of the phase operator may be evaluated 
from known properties of the creation, destruction and number operators. 
Inverting, we obtain 




cfi = exp(-/<D)(fA£ + l) 



[III-5] 



Susskind, L. and Glogower, J., Physics, 1, 49 (1964) 
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exp(/4>)= (fAf+i) y a 

/ +/ n r A~X [III-6] 

exp(-/$)=a T (iA( J + l) 

and since a (ft = fA£ + 1, it follows that 

exp (i 4>) exp( -i 4>) = 1 [ HI-7 ] 

but only in this order! Operating on number states with the phase operators, 
we obtain from Equation [ 1-26 ] 

exp(/^)|n}= (fA£ +1)" /2 a |n)=(fA£ + l)~ /2 (n) /2 |n-l)= |n-l) 

f , f „ [in-8] 

exj(-/$)|n)= a T (fA£ + l)~ /2 |n)= a T (n + l)~ /2 |n) = |n + l) 

Consequently, the only nonvanishing matrix elements of the phase operator 
are 

{n -l|exp(z <£) |n) = 1 

[ HI-9 ] 

[n + l|exp(-z <t>)\n} = 1 

The phase operators defined by Equation [ 111-36 ] do have the felicitous or 
classically analogous property of revealing magnitude independent 

information, but unfortunately they are nonHermitian operators — i.e. 

{n - 1| exp(z <£) \n) * (n |exp(z <£) \n - 1)* 

-- and, hence, cannot represent observables. However, they may be paired 

into operators that are observables - viz. 
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coscj> = I{exp(z'$) + exp(-/4>)} 
sin$ = i-{exp(/ <£) - exp(-z <£)} 



[III- 10] 



which have the following nonvanishing matrix elements: 
[n -l|cos$|n) = (n|cos$|n -1) = i 

[n -l|sin$|«) = -(n|sin$|n-l) = — 

2 i 



[III-ll] 



These nearly commuting operators 13 may be adopted as the quantum mechanical 
operators which represent (as we will demonstrate anon) the observable phase 
properties of the electromagnetic field. 

For the Fock state: 

(n| cos |n) = (n\ sin4> \n) = [ III-12a ] 

Acos$ = Asin$= J{<«|cos 2 <$>\n) -{n\ cos$|«) 2 } = Jj^ [III-12b] 

Acos* Asm* = 1^ [III-12c] 

c. The coordinate or Schrodinger representation of state: 

Recall from Equations [ I- 10a ] and [1-31] that 



13 Also, it may be easily established that the matrix elements of their commutator are given by 



(rc|[cos$,sin$K) = yS nn ,S„ 
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Jl „ ( \m CO ^ m CO 2 ~l 

^^"•(JtTTm'J [m - 131 

Therefore, the probability P(q) of eigenvalues q for a given Fock state |«) is give 
by 

Ul \ 1/ I VI 1 L/2 I / mC0 1 r m(0 2 ~\ 1/11 

p(«)=N»)|— ^ ». [J— «J«p l -— « j tni-14] 

d. Approximate "localization" of a photon: 14 

Of course a plane wave is distributed or "de-localized" in both time and space. 
Defining the "wave function for a photon" is a task fraught with danger, 15 but the 
simpler task of defining a wave function approximately localized at a given instant 
is relatively straight forward ~ viz. 

H - - l 2 l 
1 k- k 

2. Photon States of Well-defined Phase : 

Consider the state defined by 

s 

\(p) = \im(s+iy A YexpUn <p]|/z> [111-15 ] 



14 See Section 10.4.2 in Leonard Mandel and Emil Wolf, Optical Coherence and Quantum Optics, Cambridge 
Press (1995), ISBN 0-521-4171 12. 

15 See Section 1.5.4 in Marian O. Scully and M. Suhail Zubairy, Quantum Optics, Cambridge Press (1997), 
ISBN 0-521-43458. 
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Clearly, (cp |cp) = 1 given the orthonormal properties of the number states. Essential 
question: Is this state an eigenstate of the phase operators? To answer the question we 
need to consider the following potential eigenvalue equation: 

cos<£|cp) = - lim^ + l) Yl \ ^exp[z n cp] exp[z'$] |«)+^exp[z'n cp] exp[-z <£] \n) \ [III-16a] 

I n=0 «=0 J 

Using Equations [ III- 10 ] and [ III- 10 ], we obtain 



4>|cp) =- lim(s +1) Yl j Yexp[in9]|n-l)-l-Vexp[zHcp]|H-l-i)[ 

= - lim (5 +l)~ /2 j exp(zcp)^exp[z v (p]|v) + exp[-z'(p] ^exp[z'v cp] |v) ' 



2 L ^ Y7 ^ n . TJ1 . / "L *TJ TJ.Wj [m .16b] 

= coscp |cp) 

+ - lim(s +l)" y2 |exp[z scp]|s + l)-exp[z (5 + 1)9] |s}-exp[-z cp] |0)| 

so that the state |cp) fails to be a strict eigenket of cos$ by terms that diminish faster 
than (s + 1) Yl as s — >°° . Similarly, we can see that diagonal matrix elements of cos$ 
and sin$ are given by 

(cp |cos <3> |cp) = coscp jl- lim (s + 1)~'| => coscp [III-17a] 
(cp I sin <3> I cp) =sincp \ 1- lim (s + 1)~' | => sincp [III-17b] 
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Reprise of Characteristics and Properties of Phase States: 

a. The expectation value of the number operator and the fractional 
uncertainty associated with a state of well-defined phase: 

s(s +1) 



(cp| fA^ |cp) = lim (^+ l)" 1 £n=lim + 



= lim — 

2 



fractional 
uncertainty 



{<cp|^ 2 |9)-<9K |9) 2 } 



WW |q>) 



r 

^limfs+l) -1 Vn 2 - 



lim (5+I) V« 



n=0 



lim(^+l) 



n=0 



I lim !~-(2^ 2 + 5)- — s 2 ,1 
[ S ^~L6 V ; 4 JJ 1 



lim - 



73 



[III- 18a] 



[III- 18b] 



b. Expectation value of the fields associated with a single mode: 

From Equation [ III-3a ] 



{(p|E|cp) =-2 I esin(k r-O) t+ cp) lim {s+ 1) '^(n+l) 

2 £ V 1 * „ = 

=> diverges as Js for large s ! 



[Ill- 19] 
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c. Phase of field associated with single mode: 



(cp | cos <3> I cp ) = coscp 
(cp| sin$ |cp) = sincp 



[ m-20a ] 



Acos<£ = Asin$= 



<cp| cos 2 <3>|cp)-<cp|cos«3> |cp> 2 } =0 [ III-20b ] 



d. Probability of photon number: 

Finally, we may easily deduce the probability of finding n photons (i.e. the photon 
statistics) in a particular state of well defined phase — viz. 



We see that there is a equal, but small probability of any number: this agrees with the 
intuition that the magnitude of the field is completely undetermined if the phase is 
precisely known! 

3. Coherent Photon States: 16 

It would, indeed, be useful to have eigenstates of the destruction operator (electric 
or magnetic field) — viz. 



P n =ln\<?f= lira (s + l) 



[ 111-50 ] 



a i\ a i/= a i\ a 



\ 

k/ 



[111-51 ] 



Reprise of Characteristics and Properties of Coherent States: 



a. The Fock state representation of the coherent state: 



The coherent state is a Harvard invention! See R. J. Glauber, Phys. Rev. 131, 2766 (1963). 
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Since, a ^\n)= Jn+1 \n + l) and a a ^ = 9£ +1 , then (n\ a = Jn + l (n + l\ and we 
are able to write a representative of the sought state in the number state basis — viz. 

{n\a \ a)= Jn+l (n+l\a) = a(n\a) [III-52a] 

a , , . a" 



or (n|a)=-y=<n-l|a)=-p(0|a) 
Using the expansion of the identity operator, the eigenket becomes 

|a> = E|n>(n|a> = <0|a)E-| T |n>. 

n n V • 



To normalize the eigenket write 

(a|a) = (a|0}<0|a)E^-^ 



(a|0)(0|a)exp[|a| 2 ]. 



r i 



Finally, we see that 



so that { cc| 0) = <0| a) = exp \- - |a| 

. . r 1 , ,2~l V a " i \ 
|a}=exrj--|a| j 2^-7^1 I") 

n ** 



[ III-52b ] 



[ 111-53 ] 



[ 111-54 ] 



[ 111-55 ] 



is a normalize representation of the eigenkets of the destruction operator. 
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b. The expectation value of the number operator and the fractional 
uncertainty associated with a coherent state: 

(cc^Vja) = M 2 [III-56a] 



[uncertainty] <ct|5\£ |cc) |aNl n M )U n! n 



I I 2k 
a 



= ^ exp(- N 2 )ltf[n(n-l) + n]-|ar 



[ III-56b ] 



=iar 



Thus, we see that the fractional uncertainty diminishes with mean photon number! 

c. Expectation value of the electric field associated with a single mode: 

From Equation [ III-3a ] 



(cc|E|cc) = -2 / — — e|a|sin(k-r-G) t+&) [IH-57a] 

where oc= |a| exp(z't}) . 

AE = J{(a|E.E|a)-(a|E|a) 2 } = " [ ni-57b ] 



17 Similarly AH= — ^— I - — for the coherent state, so that AE AH = C fl 0)/2 V . 

c u. V 2 e V 
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d. Probability of photon number: 

From the representation of the coherent state given in Equation [ 111-55 ] we may 
easily deduce the probability of finding n photons (the photon statistics) in a 

particular coherent state is given by a Poisson distribution characterized by the 

- i i 2 

mean value n = |a| . — viz. 



P„ = ^|c^ 2 = exp[-|af 



|a| 



In 



[ 111-58 ] 



Sample Poisson Distributions - Coherent State Photon Statistics 



0-4t r\. n = kr = l 




0.2t 



J ~ r 


i = |«f = 1 






n = |ap = io 






\ -,r 














n = |ctp = i 






-i 



40. 



SO. 



120. 
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e. Phase of field associated with single mode: 



(a|cos$|a) = ^ exp 



1 

2 



1 ,..,0 



a a" + a a 



= |a|cos$ exp 



2h 



[ III-59a] 



Unfortunately, it is not possible to evaluate this summation analytically. However, 
Carruthers 18 has given an asymptotic expansion which is valid for a large mean 
number of photons ~ viz. 



{a | cos | a) =cosi3 



1 



1 



8|u| 



7+- 



|a| »1 



[ III-59b] 



f. Coherent states as a basis: 

As we will see presently, the coherent states are very useful in describing the 
quantized electromagnetic field, but, alas, there is a complication — the coherent 
states are not truly orthogonal! From Equation [ III-6 ] we see that 



(P|a) = exp[-I|af-I| P |^ ^ 
= exp[-i|af-i||3| 2 + ap*] 



[ 111-60 ] 



so that 



18 Carruthers, P. and Nieto, M. M., Phys. Rev. Lett. 14, 387 (1965) 
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(ccp) (P|cc) = exp(-H 2 - 1(3| 2 + a(3* +a* pj 

= exp(-(a-(3)(a*-|3*))=ex I (-|a-(3| 2 ) 



[ 111-61 ] 



That is, the eigenkets are approximately orthogonal only when |a-(3| is 
large ! 

g. The "displacement operator:" 

There are a growing and significant set of applications where it is useful to express 
the coherent states directly in terms of the vacuum state |0) . If we use the number 
state generating rule 



\n) = 



a 



f 



|0> 



— i. e. Equation [ 1-27 ] — the coherent state may be written in the form 

( 1 2^ [aa-) f . i 2 \ 

|a}=exp--|af I> f- 10)= exp a a T - - |cx| 2 |0) [111-62] 

V 2 J n n \ \ 2 J 



If we make us of the Baker-Hausdorff theorem, 19 we may easily show that 



The Baker-Hausdorff theorem or identity may be stated as 

exp[^ }= expjj?} exp{$} exp{-^[j? ,<B ]} 

when , [ A ]] = [$ , \A , $ ]] = 0. For a proof, see, for example, Charles P. Slichter's Principles of 

Magnetic Resonance, Appendix A or William Louisell's Radiation and Noise in Quantum Electronics. 
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|cc) = A t (a) |0) = exp(a a - a* a )|0) 



[ 111-63 ] 



4. 

SO that ^ 1 (a) may be interpreted as a creation operator which generates a 

4. 

coherent state from the vacuum. (Its adjoint operator ft (oc)= A (-a) is a 

destruction operator which destroys a state). In some treatments !A. (a) is 
described as the "displacement operator" (written £> (oc) ) 20 and the coherent states 
are called the "displaced states of the vacuum." 21 



To explore this point of view (and to give some meaning to the phase of the 
coherent state eigenvalue), we may express |a) in a two-dimensional, 

dimensionless "phase space" representation. To that end, following Equation 
[ 1-16 ], we write the dimensionless coordinate as 



= 



h 



q = a ^ exp [i y] + a exp [- i y] 



[ III-64a ] 



and the dimensionless momentum as 



( o V s 



% = 



\tn h(0 j 



■p = a ^ exp [ i ( y+ 71/ 2)] + a exp [- i ( y+ tc /2)] [ III-64b ] 



so that 



[0,7t]= 2i 



a , a 



= 2i 



[ III-64c ] 



We can (or rather you will) show that £) ' ( a ) a T> (a) = a + a and T> ' (a) a ' T> (a) = a ' + a 

See Elements of Quantum Optics, Pierre Meystre and Murray Sargent III, Spinger-Verlag (1991), ISBN 0-387- 
54190-X. 
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and since these variables are canonical 22 

((AG) 2 ){(Att) 2 }>1 [III-64d] 

Since 



f 1 

a ' = - (9-1 7l) exp[-z y] 

a = - (9 +i 7l) exp[z y] 
2 



[ 111-65 ] 



the mode field (see Equation [II- 24a]) b 

E(r,z)=z'e£ja exp[z'k-r- /to?] -a ^ exp[-z kr+z'co ?]j [ IH-66a ] 



becomes 



E(r,?)=-e e{ti cos(k r-to ?+y) + 6 sin(kr-co ?+y)} [ III_66b ] 

Since p has a coordinate space representation -i hd/dq = -i(h(ti/2) 2 d/dQ 
and q has a momentum representation ihd/d p = i (h/2 co) /2 d/d%, 23 

a a ^ - a* a -a r [a ^ - a i a,, [a ^ + a j 

=-[a r d/de+a t d/dn] [ III-67a ] 



22 



Of course, in general {(AA)' )((AB) 2 )>"^[ A,B ] )f where ( ( A A)' )=(A 2 )-<A )' 



23 If this unfamiliar, see Equations [ 1-20 ] and [ 1-22 ] in the lecture notes entitled The Interaction of Radiation 
and Matter: Semiclassical Theory. 
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and 



A ^ (a) = exp^a a ^ - a* a exp[-(a r d/dQ + a i d/dn)] [ Hl-67b ] 

Thus, $C (a) defines or generates a two-dimensional Taylor expansion when it 

acts on a function of 9 and %. In particular, if we take the "phase space" 
representation of the ground or vacuum state |0ji) as the product of two Gaussians 

f 

(see Equations [ I-10a ] and [ 1-29 ]), then A (a)|07t) represents a shift or 
displacement of this "phase space" representation - i.e. 

(07t|a) = (07t|^ ^(a)|0)= w c (e-a r ) u G (%-a) [111-68 ] 

In light of Equation [ II-23b ], |oc (?))= |cc exp(-z' (0 ?)) we can write 

[<d%\a{t)) = w c (e-|a|cos(co t+§)) t/ c (ji:- |cx| cos(co f+<|))) [111-69] 

where a= | a | exp(z (])) . 

h. The diagonal coherent-state representation of the density operator 
(Glauber-Sudarshan P-representation) : 



It may be easily established that 
= 1 



1 = ~ If IP) (PI d 2 P= jj |P) (P| d Re((3) d Im((3) [ 111-70 ] 



so that it seems quite reseasonable to look for a representation of the density matrix 
is the form 

p= lk(P)lPXPM 2 P [ni-71] 

For a pure coherent state, P is clearly a two-dimensional delta function 
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Example 1 -- Coherent state 
P ( (3 ) = 5 (2) ( (3- a ) =8 (1) ( Re ((3) - Re (a) ) 8 (1) ( Im ((3) - Im(cc) ) [ IH-72 ] 

In general, using Equation [ III-60 ] — i.e. 

r i 2 i 2 *1 

((3|a) = exp^- - |a| - - p| + a p j [ ni-60 ] 

we may find a simple procedure for finding the P-representation by writing 
{-a|p|a)=JJp((3)(-a||3)((3|a}J 2 |3 

~ fill 73 ] 

= exp(-|a| 2 ) JJ[p( P) exp(-p| 2 )jexp[a(3*-(3 a*] d 2 (3 

Thus, (-a | p | a) exp f— |a| j is the two-dimensional Fourier transform Of the 
function P( (3 ) exp ( - f3| 2 ) and we may write 

p( p )= — exp( |pf) JJ (-a| p | a) exp( |a| 2 ) exp[-a|3*+|3 a*] J 2 a [ HI-74 ] 

As a second example, consider a thermal radiation field described by a canonical 
ensemble 



cxp{-^/k B T) 
Tr[exp(-^M fl r)] 



[111-75 ] 



where 9{ = hto\a^ a + ^\- Thus, 



p=E|i- 



exp 



ftco S\ 



exp 



WW 



[111-76 ] 
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and 



[n) = Tr^p a a j= 2^ ex P 



h(0 



\ 1 
1 



[111-77 ] 



so that 



/ n \" 



n (l + (n)) 



n) In 



Thus, we can write (n| p | n) = 



and 



(!+{»»' 



6i>" 



(-a| p |a)=£ y (-a|n)(n|a) 
- (! + <«)) 

exp (-|a| 2 ) (-|a| 2 ) T ( n ) T 



exp 



!+<„) 



exp 



-|a| 



[111-78 ] 



[111-79 ] 



[111-80 ] 



Finally, we see that 

Example 2 -- Thermal radiation - a chaotic state 



( N 2 ) 



r 



7C 2 (l + (n)) [ 

= r7^ ex p(-M 2 /{")) 



1 



\1 



exp(-(3a* + a|3*) J 2 13 



[111-81 ] 



As a third example, consider Fock or number state. From Equation [ 111-55 ] we 
see that 
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exp(-|a| 2 ) „ 
{-a|p [a) = (-a[n)(n|a)= -|a| [III-82a] 



and 



P(P ) = 1 A ex P ( |(3| 2 ) jj(-H 2 )" exp[-a(3*+P a*] J 2 a 

/ ||2 \ [III-82b] 

so that 

Example 3 -- Pure Fock or number state 

exp( B| 2 ] -)2« 
^(P)=^7^-|--;5 (2) ((3) [III-82b] 
ni 3(3 3p 



i. The Glauber-Sudarshan-Klauder "optical equivalence" theorem: 

Suppose we have some "normally ordered" function 



F (N) 



f 



a , a 



v 



r VEE 



[111-83 ] 



The expectation value is given by 



Pf 



(N) 



a , a 



v 



[111-84 ] 
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Using Equation [ III-7 1 ] we see that 



/f(N) 



a , a 



1 



Tr 



n m 

= jj^(«)EEc„ m (a|a , "a">)^ 2 « 

n m 

= l\p{<*)H,c nm a"a"'d 2 a 

n m 

or, finally, the "optical equivalence" theorem 

(f (N {a,a t )=JJp(a)f (N) (a,a) 



j. The Uncertainty Relationship for {0,7t} 



Since 



a , a 



f 



: 1 we see from Equation [ III-64a ] that 



r° / \° / vv 



|a ^ a exp[2z y]+(a a) exp[-2 z yj + ^a ^ a^+^a a 



[III-85a ] 



[III-85b ] 



-la ^\ exp[2z y]-(a) 2 exp[2z y]-2^a ^\(a) 



[ 111-86 ] 



= (:A0 2 :) + 1 



where (: A :) symbollizes the normally ordered expectation value of the operator 
A . From Equation [III-85b ] 

(:A0 2 a,a^ = JJ P(a) A0 2 (oc,a*)d 2 a [111-87] 



V 
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(\ AG 2 la,a^= JJp(a) [A a* exp( / y) + Aa exp(-/ y)] 2 d 2 a [111-88 ] 



7 



/ f 

If we choose y (and P( a) ) such that ^:A0 2 



then (A0 2 )>1 and 



( A % 2 }> 1 (squeezed states)! 
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